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SUMMARY 


Studies of fission-product gamma emanation shortly after fission (1 to 
100 sec) have been initiated. The half-lives that will be established in 
this program, when used in conjunction with initial yields, will provide 
the fission-product yield history for calculational models. The studies 
undertaken include (1) an investigation of the transitive behavior of fission 
products from approximately 1 sec after fission, (2) experiments to 
establish recoil range as a parameter to use for identifying the sources of 
specific gamma rays, and (3) design of a Bteady-state tape system for 
catching and subsequently investigating recoil-fission products. This 
latter study is the result of earlier investigations and shows promise of 
providing extensive, accurate data on fission-product half-lives. 

Studies of diffusion of various nuclideB in the-Ca0-Al203-Si02 system, the 
Na 20 -Al 203 -SiC >2 system, and Small Boy soil from the Nevada test site 
have been made. The use of fission products from slow U-235 fission 
in situ has been demonstrated as a surveying method for studying diffusion 
in various matrices doped with U-235. The results of the studies are 
consistent with the compensation law, particularly with the compensated 
data developed previously in this program. The liquid-structure factor 
appears to be important in establishing diffusivities. 

Transpiration studies of the volatility of Cs from CaO-Al 203 -SiC >2 solu¬ 
tions indicate that simple gas-phase transport laws apply to transpired 
Cs. The condensed-state and vapor-Btate chemistry of this system have 
not been totally resolved. 

A mathematical description of Borption of gaseous components by a parti¬ 
cle, using condensation coefficients and condensed-state diffusion coeffi¬ 
cients, is presented. This treatment suggests that condensation coeffi¬ 
cients become important to the process when they are lower than approxi¬ 
mately 2 x 10 "k. 

Mass spectrometric studies have been made that provide estimates of the 
stabilities of the gaseous species Sb^^,, AS 4 O 6 , SbO, TeO, and Te02> 

The gaseous species As02 and Sb02 were not observed, thus their 
stabilities could not be measured. 
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INTRODUCTION 


In the past, fallout studies at Gulf General Atomic have been c oncerned 
with evaluating chemical fractionation of fission products during fallout 
formation. The studies presented in this report are an integral part of 
this evaluation, but the area to be studied has been somewhat changed. 
Henry's law studies are not, as in former studies, the only major con¬ 
cern of the work presented here. In place of Henry's law studies, mure 
effort has been spent on the delineation of the history of isotopic species 
present in the cloud. This is manifest in a program to determine half- 
lives of early fission-product nuclides. Although only preliminary studies 
have been made the results are very promising and it is anticipated that 
considerable knowledge can be obtained concerning fission-product half- 
lives in the time region of 1 to 100 sec after fission. 

More concern is being given the condensation process in fallout formation. 
Calculational studies of controlling factors for sorption of fission products 
by fallout have been made. A range of importance of condensation 
coefficients is suggested herein, and experimental measurements have 
been proposed. 

The diffusion studies program is considered to be so important that it has 
been expanded. More work on surveying methods and theoretical studies 
is presented. Studies of various matrices are important, and the matrix 
problem is of major concern in supplying data for a diffusion-controlled 
fission-product sorption model. 

Thermodynamic studies of various gaseous systems have been continued. 
These mass spectrometric studies are valuable in that they provide the 
basic data required for estimating Henry's law constants. 
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SHORT-LIVED FISSION-PRODUCT STUDIES 


As a part of this study an investigation of fission-product gamma spectra 
shortly after fission, as a function of time after fission, has been initiated. 
The purpose of this study is to delineate half-lives of some of the fission 
products that have not yet been studied in detail. This information is 
important to a model describing fallout formation because the chemical 
properties and abundances of the isotopes during fallout formation define 
the fractionation of fission products in fallout. Given initial fission- 
product yields, fission-product half-lives are used to define the quantities 
of each isotope present after detonation. 

INITIAL TRANSIENT STUDIES 

Half-life studies may be performed in several ways. A first investigation 
using a pulsing TRIGA reactor and a "rabbit" tube has been made. The 
rabbit tube is a pneumatic tube in which a sample can be rapidly delivered 
to the spectrometer-detector system from the reactor. The detector 
system consists of a 4096-channel analyzer and appropriate computational 
equipment. The pulsing-reactor/rabbit-tube system had been previously 
Bet up at the TRIGA facility. 

In this experiment, I mg of natural U in 0.2 g of pure C was subjected to 
approximately 2 x 10^ thermal neutrons/crn^ in approximately 0.02 sec, 
causing about 2 x 10^ fissions. This sample was rabbited to a position 
about one foot from a 3 cm^ Ge-Li detector and counted from the time of 
arrival (approximately 0.6 sec after pulsing) until 2 sec after pulsing. It 
was recounted during succeeding 6-sec intervals or multiples thereof for 
about 18 min, shortening the sample detector distance when appropriate. 
The resulting spectra were subjected to a computer search program to 
identify the various gamma rays according to energy and to evaluate their 
intensity. Table 1 is a listing of energies of gamma rays recovered by 
the programs during this study. To be included in this table, a peak must 
have been recovered from at least three consecutive spectra. 

None of these peaks have been identified so far. However, gamma rays 
from Cs-138 were observed in the last two spectra. Some other peaks 
have energies corresponding to reported values from certain fission 
products, but the effort afforded to this work so far has led to few identi¬ 
fications. This identification problem was previously recognized by 
Gordon, Haivey, and Nakahara.(0 
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TABLE 1 


ENERGIES OF RECOVERED GAMMA RAYS 


Energy 

(KcV) 

No. of Spectra 
Including This Peak 

Time of First Recovery 
(sec) 

85 

18 

0.6 

104 

7 

31 

Ill 

7 

258. 5 

127 

5 

381. 5 

132 

10 

0 . 6 

146 

10 

2. 5 

160 

11 

101 

169 

4 

0 . 6 

178 

6 

320 

200 

10 

101 

206 

6 

2. 5 

227 

15 

2. 5 

239 

10 

75. 5 

266 

8 

209 

292 

4 

320 

310 

6 

46 

318 

11 

101 

349 

7 

132. 5 

403 

6 

10 

412 

10 

75. 5 

438 

7 

258. 5 

458 

7 

75. 5 

464 

5 

381. 5 

544 

7 

2. 5 

590 

13 

31 

L 642 

7 

132. 5 

| 658 

3 

664. 5 

721 

8 

164 

789 

4 

443 

830 

5 

31 

912 

6 

258. 5 

968 

9 

2. 5 j 

1195 

6 

207. 5 

1236 

7 

207.5 

Y- 

3 

... " ' n ' .. ' i 

1 

/ 
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One uf the major problems with the -4096-channeL-analy/.ed data is the 
loi-.o signal-processing time (approximately 100 microseconds). This 
means that a maximum of only 10,000 pulses/sec can be received, which 

is a weak signal when divided into the individual gamma peakB. More 
extensive data will be required for accurate studies, presumably from 
steady-state system studies. Ln addition to providing more accuracy on 
the presently recovered peaks of Table 1, many more peaks will be 
recovered using a steady-state system. 

Another problem the experiment suffers from is considerable background 
signal. Bremsstrahlung is believed to be a major cause of this problem. 

RECOIL-RANGE STUDIES 


In studying fission products shortly after fission, we have encountered the 
problem of identifying the elemental sources of the gamma rays. The 
normal way of attacking the species identification problem is to perform 
chemistry to separate fission products and to then observe the decay of 
the identified, separated elemental species. Project goals of obtaining 
defined spectra between 1 sec and a few minutes after fission would seem 
to preclude, at Least for the shorter times, the chemical approach. Thus, 
an investigation oi the possibility of using the physical separation associ¬ 
ated with fission recoil has been made. 

In the first experiment, a thin fission source (dissolved normal U) cast in 
cellulose acetate (approximately 3 microns thick) backed on a 6 -micron 
Mylar film was used as a source of recoiling fission products. This film 
was sandwiched between other 6-micron Mylar films and irradiated in a GGA 
TRIGA reactor. About 1/2 hr after irradiation, or as soon as the films 
could be handled easily, the film package was undone and the two fiLms, 
one, two, and three layers ren. ved from the source film, were counted 
using the Ge-Li detector at appropriate time intervals. The source and 
film arrangement are shown in Fig. 1, and the results of these counts are 
presented as the ratio of counts in the second and third films corrected 
for radioactive decay losses (or gains) as a function of mass number. 

Before describing the results, it is helpful to understand what might be 
expected to occur according to a slightly oversimplified model of this 
experiment. First, assuming the energy of a fission-recoil product is 
deposited linearly with distance through an absorber, a particular fission 
product originating at a point source would be stopped on the surface of an 
imaginary sphere in the absorber. The radius of the sphere is the recoil 
range for the particular fission product. If this sphere is intersected by 
equally spaced planes, the areas of the portions of the recuil sphere be¬ 
tween these planes can be compared (i. e. , area is proportional to the 
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recoiled fission-product concentration). The area of the sphere between 
two planes will be the same as between another pair of equally spaced 
planes intersecting the sphere. The area between a plane inter seeling the 
sphere and one that does not will be proportional to the distance between 
the last intersecting plane and a plane tangent to the sphere parallel to 
the intersecting plane. Of course, a source plane is a multitude of source 
points with the above properties so that if the intersecting planes are 
parallel to the source plane, the above statements describe fission-product 
concentrations for this model. Finite thickness of the source plane alters 
the picture predictably. However, the plane source description is nearly 
true. In this model, then, fission-product recoil range may be distin¬ 
guished by comparing the amount of a fission product in the last film that 
catches the fission product to the amount in a film nearer the source 
plane. Fission-product recoil-energy dissipation is not truly linear with 
distance; thus, precise range determinations are more complicated, but 
the real case parallels the above presentation. 

Fission-product recoil range is mainly dependent on mass number of a 
fission product, its initial energy, and stopping properties of the medium 
in which recoiling takes place. In the nuclear scission process that 
produces a pair of fission products, each of the highly repelling nuclei 
attains a momentum of equal magnitude but opposite direction. The 
energy of fission is divided mainly into kinetic energy of these two parti¬ 
cles. Generally, the lighter-fission products made in fission are more 
highly energetic and more penetrating. While this may not be universally 
correct, it is a good first approximation. 

We have observed recoil-range variation with mass in our experimental 
measurements, as shown in Fig. 1. We believe that this curve is . 
descriptive of the particular fission-product mass and that with better 
statistics we can achieve identification of a fission-product mass from its 
experimental recoil range within one or two mass units for a considerable 
range of masses. Gamma intensity data 10 hr after fission, presented in 
this figure, fall on the curve where the statistics are good. In some 
instances this relationship appears to fail. Where it does fail (as with 
Xe-135), we have found such experimental problems as unresolved gamma 
rays, which arise from two or more fission products of significantly 
different masses. This suggests that improvements in data processing 
are required. 

A more extensive study of recoil range has been made. A series of Ge-LA 
gamma-ray spectra, made 42 min to 30 hr after a 60-min exposure of a 
sample in a GGA TRIGA reactor, has been produced. These spectra were 
obtained from recoil-fission products captured by stacked 6-micron Mylar 
films — in particular, in the second- and third-removed Mylar films from 
the U sample. In Fig. 2, a second-film-removed spectrum taken between 
26 and 66 min after fission, is presented. The identities of most of 
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Fig. 2. Spectrum for second film taken 26 to 66 min after fission 
(1 -hr irradiation) (Sheet 1 of 8) 
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Fig. 2. Spectrum for second film taken 26 to 66 min after fission 
(1-hr irradiation) (Sheet 3 of 8) 
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Fi;'. 2. Spectrum for second film taken 26 to 66 min after fission 
(1-hr irradiation) (Sheet 4 of 8) 
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g. 2. Spectrum for second film taken 26 
(1-hr irradiation) (Sheet 5 of 8) 
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Fig. 2. Spectrum for second film taken 26 to 66 min after fission 
(1-hr irradiation) (Sheet 7 of 8) 
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the peaks in this spectrum are called out, although a few of these identifi¬ 
cations are tenuous. The identifications that have been made are in agree¬ 
ment with the anticipated recoil ranges and half-lives of the fission 
products within the accuracy permitted by tht counting data. Thu data 
from this experiment (this spectrum and nine others) are presented in 
Table 2 as a gamma energy, a peak intensity ratio in the third and second 
films, an apparent half-life, and an isotope identification where an assign¬ 
ment has been made. While there may be a few assignment errors, there 
has been a great deal of success using the film intensity ratio in helping to 
establish the assignment. The problem of proper unfolding of closely 
spaced gamma lines remains, 

As a further demonstration of the recoil-range mea urements, C 3 /C 2 
ratios for 1-135 at five different times after irradiation were measured. 
The results are presented in Fig. 3. The 1-135 peak intensities were too 
weak at the first time measured and are not shown. The data at the 
second and fifth measured times arc marginal although shown on the 
graph. The deviation bars are calculated standard deviations, o , accord¬ 
ing to counting statistics. Dropping four points because of suspected bias 
(more should be dropped), 60% of the data fall within one o of C 3 /C 2 = 0.338 
and 85% within two tr, These values should statistically be 68 % and 95%, 
respectively. The very accurate data points (2. 5% it or lees) that are 
necessary for the projected studies are in good agreement with this value. 

The best data on heavy-mass-i'ission products are presented in Fig. 4. 
These data certainly suggest that recoil range can be used for nuclide 
identifications in the reported range provided the relationship between 
the nuclide number and C 3 /C 2 i H monotonic — a line is used to describe 
the data in the figure. The data define the selected line within one mass 
unit, suggesting that with good count data a nuclide's maos number can 
be identified from its C 3 /C 2 values in this region, plus or minus one 
atomic mass unit. 

STEADY-STATE SHORT-LIVED FISSION-PRODUCT STUDIES 


The requirements for making an extensive, accurate study of short-lived 
fission products can be enumerated. First, one must observe a sufficient 
number of fission-product decay events at several times after fission so 
an accurate mathematical description of birth and decay of fission prod¬ 
ucts can be obtained. Of course, sufficient energy resolution of these 
gamma counts into a specific disintegration process must be possible. 
Second, a method of associating particular isotopic disintegrations with 
each observed set of count data associated with a particular gamma 
energy must be available. The first requirement is best satisfied by 
making steady state studies of fission-product decay using a Ge-Li 
detector and an appropriate multichannel analyzer (4096 channels). Using 
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TABLE 2 


FISSION-PRODUCT GAMMA-ENERGY/RECOIL-RANGE CORRELATION 



































TABLE 2 (continued) 


Energy 

(KeV) 

c 3 /c. 

*1/2 

(hr) 

Isotope 

Energy 

(KcVj 

CJC, 

*1/2 

(hr) 

Isotope 

496 

0 . 82 


Ru-103 

266 

0. 9 

10 

Y-93 

488 

0 . 1 

2 


253 


Long 

Zr-97 

468 

0 . 8 

4. 4 

Ru-105 

249 

0 . 26 

7 

Xe-2-135 

411 

0 . 22 

Short 

Ba-141 

240 

0. 9 

3 


452 

0 . 22 

Short 

Pr-146 

228 

0. 388 

Long 

Te-132 

447 

0.8 

3 

Y -92 

220 

0. 3 

6 . 8 

1-135 

433 

0. 3 

1. 9 

Te-1-133 

211 

0. 25 

1.3 

Nd-149 

430 

0.6 

1 . 6 

Sr-92 

201 

0 . 39 

0 . 8 

Te-134 

42 3 

0. 3 

2. 3 


190 

0. 29 

0. 5 

Ba-J 41 

417 

0. 38 

7 

1-135 

181 

0. 42 

1 


405 

0. 5 

1 . 2 


167 

0. 25 

1.6 

Ba - 139 

392 


1. 4 

1-134 

156 


2 


385 



1-132 

150 

0. 4 

0.9 

Te-1-131 

364 

0. 5 

Long 

1-131 

142 

0. 85 

Long 

Tc-99 

357 

0 . 88 

Short 


129 

0. 5 


Rh-105 

354 


13 


116 


2 


343 

0. 4 

Short 


105 

0 . 1 



318 

0.65 


Rh-2-105 

92 



Nd-147 

311 

0.5 6 

1. 3 


86 


2 


306 

0. 71 

0. 5 

Tc-101 

81 


Short 

Xe -133 

293 

0.21 

Long 

Ce-143 

76 


2 


277 

0.3 5 

1 . 0 

Te-134 

59 


30 


274 

0.4 

17 


50 

0. 3 

5 
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Fig. 4. Fission-product recoil-range parameter C^/C^ 
in Mylar vs, atomic mass 
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a steady-state system one wouLd be able to observe the gamma radiation 
from fission products a set delay time after birth for a long period. This 
observation can be accomplished by allowing recoil of fission products to 
separate fuel and fission products, catching the recoils in some medium 
and studying the decay of these recoiled fission products. The second 
requirement can be satisfied by using measured recoil ranges associated 
with the established gamma energies to aid in identifying the fission- 
product emitting the observed gamma ray. 

A system that should make measuring steady-state recoil range and half- 
lives possible is shown diagrammatically in Fig. 5. Two thin Mylar tapes 
pass in front of a recoiling fission-product source. To determine the 
recoil range associated with a particular gamma ray, the counting rates 
in the two films need to be well established at a particular time after 
fission. This will be done by passing each film containing recoiled 
fission products from the steady-state fission source in front of a detec¬ 
tor at a specific delay time after fission and continuously counting more 
or less of the moving film until statistically adequate count rates of 
fission products can be established. Once recoil range and, accordingly, 
gamma-ray identifications are made, one of the recoil target films will 
be counted at various delay times after fission to provide a measure of 
specific short-lived isotope concentration versus time. From these data, 
half-lives of the short-lived fission product (1 to 100 sec) can be derived. 
From these measured half-lives and initial yields, the quantities of 
individual fission products as a function of time will be better established 
for fractionation estimations. 

Equipment for these studies, involving either neutron-beam tube fission 
studies or Cf-252 spontaneous fission studies or both, has been designed 
and will be built and used as part of a new contract. 

UNFOLDING GAMMA-RAY SPECTRA 


The various gamma rays of fission products shortly after fission are so 
abundant that measured spectra lines exhibit considerable overlap, even 
with an energy resolution of 2. 8 KeV (energy band width at one-half peak 
intensity). Most of the data reported here had resolution of that order, 
as shown in Fig. 2, The studies made so far require complex computer 
processing for roughly one-half of the observed peaks to obtain reliable 
estimates of total counts under a peak. For the purpose of obtaining 
reliable total counts from these overlapping spectra, a portion of our 
effort has been expended to establish a least-squares fit to a gamma-ray 
spectrum in the form of sums of Gaussian distributions and an exponential 
background function. The system has provided reasonable description of 
some peak complexes, but it requires optimization and further testing. 
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The data-processing system presently used is not entirely satisfactory 
for the work we are trying to accomplish, particularly since peak areas 
accurate to better than ± 5% are essential to the recoil-range identifica- 
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DIFFUSION OF RADIONUCLIDES IN MOLTEN SILICATES 


The diffusion of fission-product nuclides in molten silicate matrices is of 
major importance in any realistic model for fallout formation. ( 2 > 3) We 
have continued our investigations of this phenomenon. Studies of the dif¬ 
fusion of radiosodium and radioantimony in molten matrices from the 
sodium disilicate-albite system have been completed. An investigation of 
the diffusion of radioantimony in molten Nevada soil has been performed. 
The simultaneous diffusion of several radionuclides in molten media from 
the CaO-Al203-Si02 system has been studied using two different tech¬ 
niques. The latter study was made possible by the availability of a 4096- 
channel gamma analyzer and a high-resolution Ge-Li detector. The 
experiments have helped to clarify the role of silicate composition in 
diffusion. 

DIFFUSION IN THE Na^-Al^-SiQ^ SYSTEM 

The use of a condensed-state diffusional model to predict fractionation in 
fallout requires knowledge of diffusion rates in fallout matrix material. 
Most of the diffusion studies performed at GGA have been concerned with 
the Ca0-Al203-Si02 system. Methods of predicting diffusivities have 
been developed and tested for this system. For estimating fractionation 
during fallout formation, it would appear necessary to at least test diffu¬ 
sion-coefficient estimation methods on other systems. This has been 
partially accomplished in studies of the Na 20 -Al£ 03 -SiO^ system. 

For these investigations, four compositions from the Na^O-AI 2 O 3 -SiC >2 
ternary system were chosen. These silicates were 100% albite 
(Na 20 • Al^Oj • 6 SiO^), 100% sodium disilicate (Na^O • 2 SiO^), and 
38:62 and 60:40 mixtures of these silicates. This particular system was 
chosen for several reasons. Albite, a common mineral, can be consid¬ 
ered representative of a medium refractory soil for fallout simulation, 
and sodium disilicate is representative of a low refractory soil. Also, 
there are data reported for the leaching of radiocesium-diffused albite. ^ 
From the latter, and using the results of the present study for Na, pre¬ 
dictions of diffusivities of several other radionuclides can be made for the 
molten albite composition. This will be discussed in detail below. 

Finally, this system was chosen to study the problem of precrystalline 
structure, or configurational entropy effect, in molten silicates. 
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The first part of this study was concerned with the self-diffusion of radio- 
sodium in Lnti molten Iiidli xcca, me »UU1 eta liiai wki e uaeti iui liic planer 
source technique(3) were prepared by activating the silicates in a GGA 
TRIGA reactor to yield Na-24. Six experiments were completed in the 
temperature range of 1072' to 166 1° K. These points are superimposed 
on the phase diagram reported by Schairer and Bowen(5) that is approxi¬ 
mately reproduced in Fig. 6. The first four experime. ts with the 38:62 
and 60:40 matrices were completed. These four di r< usivity values wer» 
then combined with the results for the diffusion of Na in sodium disilic; te 
that were reported by Malkin and Mogutnov(6) to establish the compen a- 
tion law for this system.(3) The latter was obtained by treating the data 
by the method of least-squares to yield 

log 1Q D o = - 4. 33 + 0. 123 E* , (1) 

where Do and E* are the preexponential factor (cm^/sec) and activation 
energy (kcal/mole), respectively, in the Arrhenius equation for diffusion. 
An experiment was then done with sodium disilicate to compare with the 
results of Malkin and Mogutnov. At 1248°K, the oerimental value was 
found to be 1. 4 x 10-5 cm^/sec, which compareu well with the value of 
1. 5 x 10‘5 cm^/sec calculated from the equation reported by Malkin and 
Mogutnov. An empirical correlation was then found between E* and the 
percent albite composition, and E* was predicted for the 100% molten 
albite composition by extrapolation. Using Eq. 1, D 0 was found for that 
matrix and an experiment was done at 16bl'K. The experimental value of 
the diffusivity was found to be 1.6 x 10‘5 cmVsec, which compared well 
with the predicted value of 2. 5 x 10"5 cm^/ sec. This result substantiated 
the empirical correlation mentioned above. The experimental results are 
summarized in Table 3. 


TABLE 3 


DIFFUSIVITIES FOR TRANSPORT OF Na-24 IN 
MELTS FROM THE SODIUM DISILICATE-ALBITE SYSTEM 


T 

( 8 K) 

Composition 
(% Albite) (a) 

D x 10b 

(cmV sec)(b) 

1072 

38 

2. 8 

1248 

38 

6. 9 

1365 

60 

9.8 

1248 

60 

5. 0 

1248 

0 

14. 

1661 

100 

16. 


( a ) The uncertainty is approximately 2% 

(b) The uncertainty is approximately 50% 
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The Arrhenius parameters for this system are summarized in Table 4. 


TABLE 4 

ARRHENIUS PARAMETERS FOR DIFFUSION OF Na-24 IN 
MELTS FROM THE SODIUM DISILICATE-ALBITE SYSTEM 


Composition 
{% Albite) 

* 

E 

(kcal/ mole) 

lo8 IO D o 

(D 0 in cm^/sec). 

0 

11.9 (a) 

- 2. 7 -A 1 

38 

13. 5 ± 1. 0 (b) 

- 2. 67 (d) 

60 

19. 5 ± 1. 0 (b) 

- 1.93 (d) 

100 

32. 3 ± 1. 5 (c) 

- 0. 36 (C) 


^ From Malkin and Mogutnov^ 

^ Uncertainties estimated from Eq. 1 

^ From empirical correlation and Eq. 1 

(d) ^ 

From Eq. 1 


Using the results obtained for the diffusion of Na in the molten albite 
together with the leaching data for the radiocesium-diffused melts of 
albite composition reported by Lane, (4) it was possible to predict diffu- 
sivities for several species in this matrix. Lane allowed Cs to diffuse 
into the matrix, at different temperatures above the liquidus. He then 
leached the samples using 0. 1 N HC 1. From his data E* for the diffusion 
of Cs into the matrix is estimated to be 36. 6 kcal/mole. It has been 
shown(3) that E*(or logioDo) I s a linear function of the reciprocal ionic 
radius, r~l, 

* - 1 

E - a + br (2) 

for the alkali metals and I in those cases where sufficient data exist. For 
the present case, the coefficients in Eq. 2 were determined using the 
values of E* for Cs and Na. The values of E* for I, Rb, K, and Li were 
then calculated. Equation 1 was then used to calculate the corresponding 
values of logloDo* These results are summarized in Table 5. 
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TABLE 5 


EXPERIMENT A L AND PR.EDICTED ARRHENIUS PARAMETERS 
FOR DIFFUSION IN MELTS OF ALBITE COMPOSITION 


Nuclide 

* 

E 

(kcal/molc) 

lo «10°o 

(D in cm /sec) 
o 

I 

37. 8 

0. 32 

Cs 

36. 6 

0. 17 

Rb 

35. 8 

0. 07 

K 

35, 1 

- 0. 01 

Na 

32. 3 

- 0. 36 

Li 

26. 5 

- 1. 07 


It is emphasized that the values reported in Table 5 are estimates for all 
of the nuclides except Na. The estimates are, however, based on a 
rational approach to the problem of predicting diffusivities. It is antici- 
pated that the prediction will be tested experimentally for at least one 
nuclide. 

A study of the diffusion of radioantimony in molten media from the 
sodium-disilicate-albite system has been completed. The temperature 
range for the experiments was 1123® to 1687°K. Experimental techniques 
were identical with those described for the diffusion of Na in this system. 
The results are listed in Table 6. 

The character of the soil that might be involved in a near-earth detonation 
is critical to fractionation calculations. If condensed-state diffusion is a 
limiting factor in the formation of fallout particles, and we believe it is at 
least one of the important factors, then orders of magnitude in the d if - 
fusivity as demonstrated in the data of Table 6 become important. The 
data are also presented in Fig. 7 where the diffusivity is plotted on a 
logarithmic scale as a function of the reciprocal temperature. 

The experiment at 1344'K using the eutectic matrix was performed, and 
then, using the compensation law as defined by the data for the albite and 
sodium-disilicate com,positions, the diffusivity in the eutectic matrix at 
1120®K was predicted. The predicted value was identical to the experi¬ 
mental value of 2. 1 x 10-8 cm^/aec at 11 23® K. The result was consistent 
with the assumption of compensation for this system. The results are 
summarized in Table 7. 
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Diffusion coefficients for the transport of Sb-124 in 
liquids from the sodium disilicate-albite system as 
a function of temperature 








TAB LE 6 


DIFFUSION COEFFICIENTS FOR THE TRANSPORT OF RADIOANTIMONY 
IN LIQUIDS FROM THE SODIUM DISILICATE-ALBITE SYSTEM 


, Matrix 

T 

(" K) (a) 

D x 10 8 

, 2/ , (b) 

(cm /sec) 

Albite 

mm 

55 

Albite 

XU 

25 

Albite 

1859'- 

15 

Albite 

mSSSm-: 

1. 9 

Sodium Disilicate 

1495 

410 

Sodium Disilicate 

1407 

170 

Sodium Disilicate 

1301 

37 

Sodium Disilicate 

1195 

5. 1 

Eutectic 

13 44 

86 

Eutectic 

1123 

2. 1 


( 21 ) 

The uncertainty is approximately 10°K 
' ' The uncertainty is approximately 20% 


TABLE 7 

ARRHENIUS PARAMETERS FOR Sb-124 DIFFUSION IN MELTS 
FROM THE SODIUM DISILICATE-ALBITE SYSTEM 


Composition 
(% Albite) 

* 

E 

(kcal/ mole) 

lo ® 10 D o 

2 / 

(D in cm /sec) 
o 

0 

5 2. 2 

2. 29 

38 

50. 8 

2. 20 

100 

87. 6 

5. 0b 


In Fig. 8, the energies of activation for diffusion of Na and Sb are com¬ 
pared with the sodium-disilicate phase diagram. It is emphasized that all 
of llie&e experiments were conducted in the liquidus region. For both Na 
and Sb it appears that there is at least a qualitative reflection of the phase 
diagram by the activation energy. This effect, although evident, is not as 
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pronounced as that found for the SbOx*CaO-Al2C>3-Si02 system. (3) How¬ 
ever, the data show evidence of precrystalline structure in the liquidus. 
This effect has been observed in several other silicate systems, for 
example.in the Na20-SiC>2 system by Dertev and Voronkova.(') Their 
data for the diffusion of Na in these molten silicates reflect the equilib¬ 
rium phase diagram; i. e. , maxima for the diffusivity appear at two 
eutectics and at the composition corresponding to NazO • 2 SiO^. This 
interesting problem deserves further attention and we anticipate conduct¬ 
ing experiments to further define it. 

DIFFUSION IN MOLTEN NEVADA SOIL 


An investigation of the diffusion of radioantimony in molten Nevada soil 
taken from the Small Boy detonation site was completed. The purpose of 
the study was to make a second comparison of Nevada soil with the 
Ca0-Al203-Si02 eutectic, which is the reference matrix used in the GGA 
laboratory. 

Sixty-day Sb-124 of high specific activity was purchased from ORNL. The 
chloride solution was converted to a nitrate solution by several treatments 
with nitric acid followed by evaporation to near dryness. This prepara¬ 
tion was added to powdered Nevada glass, melted, redivided, and 
remelted to ensure homogeneity. The final concentration of the redivided 
source was roughly 10-2 w t-% Sb. As in similar previous experiments,(3) 
an effort was made to obtain a minimum of four points and an order-of- 
magnitude change in the concentration gradient for each isothermal exper¬ 
iment, and a minimum of four points to define the Arrhenius coefficients 
for a particular system. 

The Nevada matrix used for diffusion measurements with radiocesium(2) 
was finely divided and heated at approximately 1400°C for several days. 
This was done in an effort to improve the vitreous quality of the matrix, 
since the first preparation contained a few small bubbles and crystals. (2) 
The glassy character of the silicate was improved by this treatment and, 
following microscopic examination, a representative sample was analyzed 
using X-ray diffraction. The diffraction pattern was compared with both 
the patterns of the original soil sample and the soil that was preheated at 
1000° C. Although a few lines were perceptible in the diffraction pattern 
of the "glass, " they were faint and diffuse. The resolution of the camera 
for this work was roughly 50A. The crystallites responsible for the lines 
could not be ascertained although a search of the ASTM index was made. 

Experiments with the diffusion of Sb-124 in the Nevada matrix were done 
in the temperature range of 1218° to 1766° K. The results are listed in 
Table 8. 
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TABLE 8 


DIFFUSION COEFFICIENTS FOR THE TRANSPORT OF 
RADIOANTIMONY IN THE NEVADA MATRIX 


T 

(°K) (a) 

D x 10 8 . 

(cm^/ sec)^^ 

1766 

69 

1695 

32 

1599 

13 

1499 

6. 0 

1416 

5. 4 

1218 

3. 4 


^ The uncertainty is approximately 10“ K 
^ The uncertainty is approximately 20% 


The data from Table 8 are also shown in Fig. 9 where the diffusivity is 
plotted on a logarithmic scale as a function of the reciprocal temperature. 
In this figure, it is apparent that the data fall into two sets. Microscopic 
examination of the samples showed definite signs of devitrification in the 
low-temperature samples although the high-temperature samples 
remained glassy. It was concluded that the low-temperature data involved 
diffusion in increasingly devitiified matrices. The annealing periods for 
the three low-temperature samples were 16 hr (1218°K), 3 hr (1416°K), 
and 1.75 hr (1499°K). These times are sufficiently long to allow devitri¬ 
fication if the temperatures are below the liquidus temperature. If this 
interpretation is correct, then the point of intersection of the two 
Arrhenius curves in Fig. 9, 1290° C, should correspond to the approxi¬ 
mate liquidus temperature of a "representative" glassy Nevada soil. 

There is some uncertainty associated with the assumption of Arrhenius 
havior for the low-temperature data, which are summarized in the fol¬ 
lowing equations: 


Sb/devitrified Nevada: 

log io D = 

- 6. 08 - (1- 70 x 10 3 /T) . 

(3) 

Sb/vitreous Nevada: 

log 10 D = 

1. 48 - (13. 5 x 10 3 /T) . 

(4) 
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Equation (4) is to be compared with the equation for transport of radio¬ 
antimony in the 1450° K eutectic of the CaO-Al-,0,-SiO-, ternary system, O) 

<- J c. 


Sb/ 


vitreous eutectic; 


= in - (14. 1 X 10 /T) . 




It appears that the transport of Sb in the vitreous Nevada matrix is simi¬ 
lar to that in the eutectic matrix. This result strengthens the statement 
that the CaO-AI 2 O 3 -SiO 2 eutectic matrix is a fair model of vitreous 
Nevada soil for Small Boy calculations.^) Note that the times involved in 
cooling of an atomic-device cloud are small with respect to the time 
required to complete experiments in the laboratory, and devitrification of 
a Nevada glass would probably not occur. In this sense, extrapolation of 
the Arrhenius plot to lower temperatures from the high-temperature data 
is justified for fallout calculations if the soil has been at sufficiently high 
temperature. It is anticipated that further diffusion studies will be done 
using the Nevada matrix. 

DIFFUSION IN T HE CaO-A^O ^-SiQ^ SY STEM 

A vaporization technique has been used to study the simultaneous diffusion 
of several fission-product nuclides in the 1450“ K eutectic of the 
Ca 0 -Al 203 *Si 02 system above the liquidus. The results of these studies 
were reported in a paper submitted for publication in the Advances in 
Chemistry Series. (®) 

The diffusion of several fission-product nuclides in other silicate glasses 
from this system has been studied in the temperature range of 1508“ to 
1813°K using the diffusion-couple method. Platinum capillary crucibles 
were partly filled by melting the matrix in place. The source was then 
melted in place on top of the matrix. Sources were prepared by dissolving 
93% enriched UO^ in the matrix at the 1% level followed by activation with 
neutrons using the GGA TRIGA reactor. The sources were allowed to 
decay somewhat after activation, prior to being used. The capillaries 
were annealed at 1508° , 1683° , 1773° , and 18l3°K for known periods. 
After annealing, the capillaries were sectioned as in the plane source 
technique. (3) The capillary length was measured using micrometer 
calipers after each sectioning, and the abraded material together with the 
carborundum paper was gamma-analyzed using a Ge-Li detector and a 
4096-channel analyzer. With this analyzer, it was possible to monitor 
many individual fission-product peaks, Often more than one photopeak of 
a particular species could be studied, thus providing a measure of the 
internal consistency of the data. 
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The mathematical model used in these investigations is the solution of the 


J i---a..: ...at 


r.CCm* ai>i^u-*iiiucpcauciu uinuoivuy 


iur uiffubion 


from an infinite source into an infinite sink; i. e. , the lengths of the 
source and sink are large with respect to the diffusion distance. The 
pertinent equation is 


2 C 


C erfc 
o 


(x /2 (Dt ) 1//2 ) , 


( 6 ) 


where C is the concentration of the diffusing species, C Q is the initial 
concentration in the source, D is the diffusivity (cm^/sec) and t is the 
time. The axis of symmetry of the concentration is at the original inter¬ 
face x = 0 of the source and sink where 2C = C D . The value of D is 
obtained by plotting values of 0/C o against the corresponding values of x 
on probability graph paper and finding the slope of the best straight line. 
Experimental data conform well to the model. 

Two different matrices were used in these studies. Matrix A wt-% com¬ 
position was 41 Si 02-22 Ai 2 C >3 -37 CaO, and Matrix B was 52 Si02- 18 
AI 2 O 3 - 3 O CaO. The experiments using Matrix A were done at 1773'K. 
These results are listed in Table 9. 


TABLE 9 

DIFFUSION COEFFICIENTS FOR TRANSPORT OF 
FISSION-PRODUCT NUCLIDES IN MATRIX A AT 1773°K 


Species 

Atomic Mass 

D x 10 7 
(cm /sec) 

Te 

132 

3. 5 

I 

131 

3. 8 

Tc 

99 

3. 5 

Zr 

95 

5. 1 

Mo 

99 

3. 7 

La 

140 

6 . 1 

Ce 

143 

4. 2 


From the data in Table 9, it appears that all of the species diffuse at 
roughly the same rate with an average diffusivity of 4. 2 x 10*? cm^/sec at 
1773° K. Thus, the critical temperature T* is approximately 177 3" K for 
Matrix A if diffusion for this system is assumed to be compensated. (3) 
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Diffusion in Matrix B was studied at three different temperatures, 1313°, 
1683° , and 1508°K- The last temperature is below the liquidus tempera¬ 
ture ol this ternary composition, it was lound that by extrapolating the 
Arrhenius line, obtained from the two higher temperatures, to 1508°K, 
most of the experimental diffusivitics at this temperature were greatly 
underestimated. Some devitrification of this matrix probably occurred. 
The data at 1508°K are not considered pertinent. The data for the two 
higher temperatures are listed in Table 10. 


TABLE 10 


DIFFUSION COEFFICIENTS FOR TRANSPORT OF 
FISSION-PRODUCT NUCLIDES IN MATRIX B 


Species 

D x 10 8 
(cmV sec) 

18 13° K 

1683° K 

Xe 

9. 5 

2. 7 

Nd 

19 

4. 9 

Ce 

19 

10 

La 

30 

13 

Te 

12 

3. 1 

Ba 

79 

24 

Zr 

7. 5 

2. 5 

Ru 

28 

4. 9 


Although the data in Table 10 are preliminary, they afford approximate 
values of the coefficients in the diffusion equation, 

log 10 D = a - (b/T) . (7) 

The values of a. and b. are listed in Table 11. 
l l 

The data in Table 11 appear to be compensated, i. e. , a decreases with 
b in a fairly regular fashion. In Fig. 10, the compensation correlation is 
shown for these data. The line drawn through the points is a least- 
squares-fit to the data with unit weighting assuming a linear correlation. 
This observation indicates that the results of these preliminary experi¬ 
ments are internally consistent and that compensation prevails for 
Matrix B. It was also calculated that, on a statistical basis, the 
"critical" temperature for this system is 1808°K and the "critical" dif- 
fusivity is 1. 9 x lO"? cm^/sec, 
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Compensation correlation for the transport of 
fission-product species in Matrix B 
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TABLE 1 1 


COEFFICIENTS IN THE DIFFUSION EQUATION FOR 
TRANSPORT IN MATRIX B 


Species 

a 

b x 10‘ 3 
(° K) 

Ru 

3.29 

17. 8 

Te 

0. 86 

14. 1 

Nd 

0. 78 

13. 6 

Xe 

0. 10 

12. 9 

Ba 

0. 63 

12. 2 

Zr 

-0. 90 

11. 3 

La 

-1. 73 

8. 7 

Ce 

-3. 23 

6. 3 


Although the results reported here are preliminary, it is possible to 
compare them with previously reported results!**) for the 1450°K eutectic 
of the CaO-AI 2 O 3 -Si 02 system. In a qualitative sense, diffusion in the 
eutectic matrix is more rapid than in the more refractory matrices, A 
and B. A major purpose of this study is to attempt a quantitative evalua¬ 
tion of the effect of composition upon the simultaneous transport of 
several nuclides in silicate systems. To accomplish this task, more 
extensive data of this kind are required. 

COMPENSATION LAW 


The utility of the compensation law, or compensation effect, in predicting 
the diffusion of fission-product nuclides in silicates has been considered 
to be important in this laboratory. (<3* 3) a paper on this subject has been 
submitted for publication. (9) The experimental data obtained during this 
year using the plane-source technique!^) have been combined with pre¬ 
viously reported data(^’ 3) for diffusion in silicates obtained in this 
laboratory. The resulting set of 28 datum points has been treated using 
the compensation effect; the correlation between log[QD 0 and E >,! /4. 575, 
whore D 0 is the preexponential coefficient (cmVsec), and E* is the 
activation energy (kcal/mole) in the Arrhenius equation. The correlation 
is shown in Fig. 1 1. The curve drawn through the points is the least- 
squares-fit 

log j qD q -■ (-4. 37 + 0. 1 10 E V ) ± 0. 60 , (8) 

where the uncertainty is the standard error of estimate. The point cor¬ 
responding to diffusion of radioantimony in the devitrified Nevada matrix 
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Fig. 11. Compensation correlation for data obtained at GGA for 
diffusion in silicates. Solid circles refer to data 
obtained this contract year 
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was omitted from the least-squares treatment since it is more than four 
standard deviations from the least-squares Line, and transport can cer¬ 
tainly be expected to be mechanistically different from that in glasses. 

The importance of this correlation is that it involves the diffusion of 10 
radionuclides in glasses that range in composition Irom low-reiractory 
sodium disiiicate to high-refractory Nevada soil and CaO-A 1 3 -SiO^ The 

ranges in the coefficients are roughly 114 kcal/mole for the activation 
energy and 14 orders of magnitude for D Q . In a qualitative sense, the 
correlation can be used to roughly estimate diffusion coefficients since 
diffusion in refractory matrices and large nuclides involves large values 
of n o and E*. and diffusion in low-refractory matrices and smaLl 
nuclides involves small values. Quantitatively, the correlation can be 
used to predict diffusivities if an independent estimate of either D c or E* 
can be made. This has been discussed previously. (3) This correlation 
will be used as the reference for any further calculations done using the 
Norman^) model and the Korts and Norman! 1 ^ c.alculational scheme. 
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HENRY'S LAW CONSTANT MEASUREMENTS 


/ 2 1 1 -1 31 

In recent reports, 1 ’ } Henry's law data have been presented for 

several radionuclides in equilibrium with molten .matrices from the 
CaO - A1 2 O 3 -SiO^ ternary system. The list of radionuclides that are amen¬ 
able to experimentation with this system, using the transpiration technique, 
appears to be exhausted. The two major limitations are satisfactory 
volatility and lack of interaction with the Pt apparatus. It had been antic¬ 
ipated that a study would be performed using radiotechnetium. However, 
Tc was found to react with, or diffuse into, Pt and these experiments 
were discontinued. Further Henry's law constant measurements may be 
extended to other matrices using fission products already found amenable 
to study in Pt. 

A study was completed that evaluated the effect of gas-phase diffusion on 
transpiration experiments as done in this laboratory. As a test of the 
transpiration method used in studying Cs Henry's law constants , ^ ^ heli¬ 
um was used in place of Ar carrier-gas. This change was made for both 
the high- and low-concentration Cs samples described p» uviously. ^ 

This study was done at approximately 1200°C, Its purpose was to deter¬ 
mine if Cs was being transported to the collector in the gas phase. It is 
possible that, at the lower flow rate, transport limits might have been 
achieved by diffusion along the walls of the capillary instead of in the gas 
phase. The reported lack of a complete understanding of the 
Cs-CaO-A1 ^C^-SiO^ system,^) enhanced the importance of this 
exoeriment. 


Gas-phase diffusivities can be estimated from Hirshfelder's descrip¬ 
tion, (14) 

3/? i 

0.001858T ' [(M 1 +M 2 )/M 1 M 2 ] 2 


(9) 


where Dj ( 2 i s the interdiffusion coefficient, M is the molecular weight of 
a gaseous species, ito| 2 an arithmetic mean cross section, ilj-j is the 
collision integral evaluated at the geometric mean of kT/t x (k/( x is a 
Leonard-Jones force constant for x), T is the temperature, and P is the 
pressure. A ratio Dj j t 3 can be developed, 
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(10) 


1 , l 


D 


1.3 


J < M ! + M 2> M 3 
» (M 1 + 

, V / f /,V / 12 ' 

' T -> I I C 1 / 1 


<Ti + a 3 




where («j x /k)*/^ is taken as inversely proportional to Qgjj , a reason¬ 
able as sumption for values of kT/ «^ larger than 10 . Above 1200“C this 
assumption appears to be adequate. Cross sections and values of e x /k 
were taken from Hirshfelder, and (r^g, which is not critical, was taken 
as 4.5A. The value of Dj t z/Di , 3, where 1 represents Cs-1 37, 2 repre¬ 
sents He, and 3 represents Ar, was calculated to be 2.87. Experimen¬ 
tally comparable Dj x ratios were obtained; 2. 2 was obtained for the high 
Cs concentrations and 2.7 for the low Cs concentrations. These ratios 
are in excellent agreement with the theoretical value within the experi¬ 
mental error. Thus, gas-phase diffusion is considered to be the control¬ 
ling process in the Cs transpiration studies at low gas-flow rates. 


Those new data do not provide an explanation of the lack of oxygen or 
water-pressure sensitivities described in the last report. ^ They do, 
however, confirm the method of analysis used to describe these previous¬ 
ly reported transpiration results. 
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THE ROLE OF CONDENSATION 
COEFFICIENTS IN FALLOUT FORMATION 

To calculate the importance of condensation coefficients in fallout forma¬ 
tion, a system has been developed to evaluate diffusion absorption in 
spherical particles where a surface reaction plays a role. 

Let us assume that the surface condensation rate is proportional to the 
condensing species pressure P (atm) and reevaporation is proportional to 
its surface concentration C s (g/g). Thus, for an isothermal system, 


~dt = Kq c P - K '"e C s 1 (H) 

where C is the average concentration of the species in a particle, t is the 
time (sec), a c is the condensation coefficient, u e is the evaporation coef¬ 
ficient, and K and K 1 are constants. To simplify the system, » c and a e 
will be taken as equal to a, a constant independent of concentration or 
pressure. At equilibrium 

O • § * K«P-K'« C<J . (12) 


where C s f is the equilibrium surface concentration. Thus, since 
P = HC s f where H is the Henry's law constant (atm/g/g), 

KH = K 1 , (13) 


§ = KH “ (l- c =) 

For an initial rate of absorption (C s - 0), 



and according to the Langmuir evaporation expression 


(14) 


(15) 



44. 3 3 


(4tr r 2 ) 


P» 


\/MA 


( 16 ) 
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where m is the quantity condensed (g), r is the radius (cm) of the sphere, 
M is the molecular weight (g/mole) of the gaseous species, and T is the 
temperature {~ K). Now, 


/dm\ _ 4irr'p /dC \ 

\ dt A" 3 UA ' 

where p is the particle density (g/cm 3 ) and 


<«),■ 


132, 99 Pq. 7m/T 
P r 


17) 


18) 


anil thus, 


K = 


_ 132. 99 7m/T 


P r 


(19) 


It follows, then, that the rate of absorption of a species by the particle is 


dC _ 132. 99 Ha i/m/T /P 


dt 


pr 


«■=.) 


( 20 ) 


( 16 ) 


To relate this information to Newman's __ description of diffusion and 

surface emission (Crank ^ 17 ^), the term ^ has to be related to —, 
Fick's law for a sphere is written dr 


dC _ 
dt 


D_ d_ / 2 dC \ 

2 Or. ( r i Or. ) ’ 

r i 1 V 1 / 


[ 21 ) 


where D is the diffusion coefficient in the condensed state. Now 
3 _ r r 

4tt r dC - f a 7 9C, _ , „ f 3 / 2 9C \ 

—;— Tr “ / 4 ir r. —dr. - 4 ttD / -— |r. - ] dr 

3 dt J 0 l 8t i J Q 9r. ^ i Or. J i 


4tt Dr 


ac 

_s 

ar 


( 22 ) 
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Therefore, 
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Finally, 


( 23 ) 


ac 

S 

dr 


44. 33 Ha Vm/'T /p ^ \ , /p ^ \ 

-Dfi- U ■/ " U C s) * 


(24) 


as a description of the surface condition where condensation and evapora¬ 
tion coefficients are important. NewrnanO^) has solved the diffusion 
problem for this boundary condition and CrankO?) gives the solution as 



(25) 


where F is the fractional absorption where the sphere starts with zero (or 
uniform) concentration, and 

3 cot 6 - 1 - rh . (26) 

r n r n 


The solutions to this problem are presented by Crank for various values 
of rh. If rh is more than 10 there is little difference between the solution 
to Eqs. 25 and 26 and that for a strictly diffusion-controlled process. 
Condensation problems, then, generally become important for rh less 
than 10. Higher values of rh may become more important at low absorp¬ 
tions. However, as a semiquantitative estimate of the importance of con¬ 
densation coefficients, if H = 1, T = 1673, M = 104, p = 2. 5, D ~ 1 x 10"^, 
and r - 0. 01, then the a for which condensation problems begin (rh = 10) 
is 2 x 10~^; that is, only smaller values of a will cause condensation 
problems. The values employed were selected from the middle of the 
ranges of the quantities important for fallout calculations. However, 
these ranges are wide and, thus, condensation problems are undoubtedly 
more important than would be suggested by a c ~ Z x 10*6 limit. AdamsO^) 
has used values of a this low and lower to explain some experimental data. 


In a more exact application of Eqs, 25 and 26 to the problem of fallout, 
difficulty in calculation at low absorptions is encountered. For fallout 
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problems, important values of F are often below 0. 1. This involves 
values of n, r, D, t, and p n such that the series in Eq. Z5 is slowly con¬ 
verging. Attempts were made to transform this equation or to simplify 
it to render calculation simpler. First attempts nave been unsucceasiul. 
It appears that a computer program is required to accurately evaluate 
this question throughout the entire range of F values. 

Another planned test of this system is the use of the SLIDER Code,^ 
described at the April 1968 symposium on fallout, to calculate the effects 
of combinations of gas-phase diffusion, condensed-state diffusion, and 
condensation rate on absorption by a spherical particle. Numerical cal¬ 
culations of this kind should be useful in describing recent work by 
Adams! 18), 
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MASS SPECTROMETRY OF OXIDE SYSTEMS 


The molecular character of gaseous -fis sion products and their oxides is 
of critical importance in describing fallout formation. These considera¬ 
tions are involved in many of the mechanisms of particle-gas interaction; 
e. g. , gas-phase diffusion, condensation, Henry's law equilibria, 
condensed-state diffusion, and reevaporation. During the past year, the 
method of high-temperature Knudsen-eell mass spectrometry has been 
used to define gaseous equilibria between metals, metallic oxides, and 
O for Sb, As, and Te. 

An attempt to evaluate the stabilities of MO{g) and M02(g) (M = As, Sb) 
has been made. In particular, the reaction 

M 4 0 6 (g) rr 4 MO(g) + O z (27) 

was studied. The metal oxides were dissolved in B 2 O 3 to reduce their 
activities. These melts were placed in an Ir Knudsen cell provided with 
an oxygen inlet. The cell was heated by radiation and electron bombard¬ 
ment. Temperatures were measured by either a Pt-Pt 10% Rh thermo¬ 
couple or an optical pyrometer. The molecules M^jO^g), MO(g) and 
were the only species observed effusing from the Knudsen cell under 
these conditions. Appearance potentials of the oxides were found to be 
10. 6 , 8 . 9, 9. 9, and 9. 0 eV for As 4 0^ + , AsO + , Sb^jO^, and SbO + , 
respectively, with an estimated uncertainty of 0. 5 eV. 

Isothermal dependence of the peak intensity ratio M 4 C> 6 + :(MO + )^ on the 
0 Z + signal showed that Reaction 27 was being studied. 

Temperature-dependence measurements of the equilibrium constants for 
Reaction 27 were made. Four determinations were made for each reac¬ 
tion with an average of 7 points each. Figure 12 and 13 are typical 
temperature-dependence plots for Reaction 27 for As and Sb. 

Entropies for this reaction were calculated from equilibrium constants 
obtained from Ag pressure calibrations of the systems. Table 12 is a 
summary of the heats and entropies. The data, including the lack of 
observation of MO<>(g), suggest that MO(g) may be the important gaseous 
species for these elements during fallout formation.- 
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Fig. 1 Z. A typical temperature dependence of t.he equilibrium 
constant of Reaction LI for As 
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TAIU.E I f 


SUMMARY OF 

THE HEATS 

AND ENTROPIES 

FOR REACTION 

M 

Temp. Range 
(°K) 

AH . , 

mid -temp 

(Real/ mole) 

AS . , 

mi d • temp 

(eu) 

As 

10^5 - 1175 

5 >09 ± 4 

147 

Sb 

U85 - 1415 

3 1 ^ :fc 6 

1 28 


A study of the reaction 


SbO(g) .♦ Sb(g) + 1/iO^ 


(Z8) 


has also been completed. 

A mixture of N.i^OOj and Sb^O}, which had been heated to its melting 
point, was put into an lr Knudsen ceil. The cuLl was placed in the mass 
spectrometer and heated by electron bombardment. Temperatures were 
measured with an optical pyrometer. 

The ions 0£ + , Sb + , and SbO + were obseived after considerable CO £ was 
evolved. Appearance potentials wer«' found to be 8, 7 and 8. 1 eV for Sb + 
and SbO + , respectively. The appearance potential for Sb + is in pood 
agreement with 8. (> eV reported in Kiser's compilation. The appear 

aiuc potential for SbO + is t). 9 eV lower than that found during the study 
of Reaction i7. The value of 8. 1 eV is believed to be more reliable 
because the SbO ^ signal was a factor of 50 higher. 

Four sets of temperature-dependence measurements of the equilibrium 
constant for Reaction <f8 were made. They yielded a 

AH°i 4 i|j - 58. 1 ± 1. b kcal/mole. The temperature range of the measure¬ 
ments was 1405° to 15H/°K. A silver calibration was done and partial 
pressures were calculated. An entropy was calculated for Reaction ^8 
from the equilibrium constant obtained from the pressure calibration. 
Table 1 5 is a summary of the results. 

A AH° [ cq)o of 4t>7 kcal/mole is obtained for 

Sb O (g) -- 4 Sb(g) +3 0, (<M) 

4 o c 


SO 








TABLE 13 


RESULTS OF THE PRESSURE CALIBRATION AND 
CALCULATED ENTROPY FOR REACTION 28 


Tern]) 

P 

SbO 

P Sb 

P C2 

AS 

(*K) 

(atm) 

(atm) 

(atm) 

(eu) 

1496 

1. 13 x 10" 6 

5.63 x 10 

2. 38 x id' 6 

11.2 


by combining the results oi Reactions 27 and 28. Using Allj^oo 
60.7 kea!./mol(^) for 1/2 02 ^ 0 and 315 kcal/mole for Reaction 27 as 
reported previously, values of 98.8 kcal/mole for SbO(g) and 
831 kcal/mole for Sb406(g) are obtained for 

Gaydon^^ reports a dissociation energy, of 74 ± 9 kcal/mole at 0 U K 
for SbO(g). This value is an extrapolated spectroscopic value in which 
the ground state of the molecule was uncertain. Brewerreports an 
energy of dissociation to the atoms at 0°K for SbijO^g) of 8H(> ± 

15 kcal/mole, a difference of only 4% from the value found in this study. 

(2 3) 

The Te-O system is described elsewhere,. 
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ta abstract 


Studies of fission-product gamm» emanation shortly after fission (1 to 100 sec) have been initiated. The half-lives that will be 
established in this program, when used in conjunction with initial yields, will provide the fission-product yield history for 
calculational models. The studies undertaken include (1) an investigation of the transitive behavior of fission products from 
approximately 1 sec after fission, (2) experiments to establish recoil range as a parameter to use for identifying the sources of 
specific gamma rays, and (3) design of a steady-state tape system for catching and subsequently investigating recoil-fission 
products. This latter study is the result of earlier investigations and shows promise of providing extensive, accurate data on 
fission-product half-lives. 

Studies of diffusion of various nuclides in the CaO-Ai^OySiOj system, the NajO-AljOj-SiOj system, and Small Boy soil from the 
Nevada test site have been made. The use of fission products from stow U-235 fission in situ has been demonstrated as a 
surveying method for studying diffusion in various matrices doped with U-235. The results of the studies are consistent with 
the compensation law. particularly with the compensated data developed previously in this program. The liquid-structure factor 
appears to be important in establishing diffusivities. 

rranspirattor- studies of the volatility of Cs from CaO-AljOySiOj solutions indicate that simple gavphase transport laws apply 
to transpired Cs. The condensed-state and vapor-state chemistry of this system have not been totally resolved. 

A mathematical description of sorption of gaseous components by a particle, using condensation coefficients and condensed-state 
diffusion coefficients, is presented. This treatment suggests that condensation coefficients become important to the process only 
when they are lower than approximately 2 x \ 0 ' b . 


Mass spectrometric studies have been made that provide estimates of the stabilities of the gaseous species Sb^O^, As^Og, SbO, TeO 
and Tkij- The gaseous species AsOj and SbOj weie not observed, thus their stabilities could not be measured. 
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